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Electronic structure calculations using the augmented spherical wave method have been performed 
for CuRh02- For this semiconductor crystallizing in the delafossite structure, it is found that the 
valence band maximum is mainly due to the 4d t2 g orbitals of Rh 3+ . The structural characterizations 
of CuRhi_ x Mg x C>2 show a broad range of Mg 2+ substitution for Rh 3+ in this series, up to about 
12%. Measurements of the resistivity and thermopower of the doped systems show a Fermi liquid-like 
behavior for temperatures up to about 1000 K, resulting in a large weakly temperature dependent 
power factor. The thermopower is discussed both within the Boltzmann equation approach as based 
on the electronic structure calculations and the temperature independent correlation functions ratio 
approximation as based on the Kubo formalism. 



PACS numbers: 71.20.-b, 72.15.Eb, 73.90.+I 

Keywords: electronic structure, low-dimensional compounds 



I. INTRODUCTION 

The search for new thermoelectric materials in order 
to convert waste-heat into electricity has motivated nu- 
merous studies on transition-metal oxides. One advan- 
tage of these materials over several others lies in their 
ability to be used at elevated temperatures in air. This 
opens the route to target systems releasing heat at tem- 
peratures as high as 1000 K. Among the studied p-type 
thermoelectric oxides, the layered ones such as Na x Co02, 
misfit cobaltites, or, more recently, the delafossites, all 
with structures containing Cdl2-type layers, have been 
particularly investigated according to the richness of 
their physical properties. For instance, the thermo- 
electric performance of several AMO2 delafossites have 
been measured, leading to the following dimcnsionless fig- 
ures of merit ZT of 0.04 at 800K (CuCro.97Mgo.03O2); 1 - 
0.14 at HOOK (CuFe .99Nio.oi0 2 )r and 0.15 at 1000K 
(CuRh .9oMgo.io02)£ Their crystal structures can all be 
described as a delafossite-type, in which M0 2 dayers of 
edge sharing MC>6 octahedra alternate along the c-axis 
with layers of monovalent Cu + cations, the latter ex- 
hibiting a dumbbell O-Cu-0 coordination. 

Although the thermoelectric properties of the Cu- 
based delafossites have been measured, their origin re- 
mains a subject of controversy. In the CuCri_ x Mg x 02 
system, the measurements have been interpreted by con- 
sidering different active layers for the electrical transport 
- cither in the Cu or Cr0 2 layers - whereas in a recent re- 
port on CuRh02 , the electrical conductivities of both Cu 
and Rh02 layers have been proposed to be comparable 
at 300 

In order to shed light on the respective role of 
the layers on the transport properties in delafos- 
sites, electronic structure calculations have proven to 
be useful^LM.iSJiiSiia In the AC0O2 delafossite 



with A = Pt or Pd, these calculations have demon- 
strated that their low resistivities (~ 5 [iQ cm at room 
temperatur e) 14 ' 15 ) come almost exclusively from the in- 
plane d orbitals of the A + cationsj^i^ For CuMC>2, on the 
opposite, it is found that the ti g states of the M = Cr 3+ 
cations provide the most important contribution at the 
valence band maximum, with spin polarization supported 
by the experimental evidences for negative magnctorcsis- 
tance and magnetothcrmopoweri 1 ^ However, as has been 
recently reported for CUYO2 r where the 4c?° stable elec- 
tronic configuration of Y 3+ precludes any participation 
of Y to the transport, there exists some cases where the 
copper cations are contributing to the charge derealiza- 
tion. However, this delafossite belongs to those char- 
acterized by a large M cation (r Y 3+ = 0.090 nm versus 
r Cr 3 + = 0.0615 nm) favoring the incorporation of extra 
O 2 " anions into the copper layer so that in that case the 
physics might be different. In that respect, the physics 
of the oxygen stoichiometric CuMC>2 subclass of delafos- 
sites is of interest, especially if we consider the possibility 
to control the magnetism at the M-site. This is outlined 
by the multiferroic behavior exhibited by the CUMO2 de- 
lafossites for M = Cr 3+ or Fe 3+ , which are both showing 
electric polarization induced by incommensurate antifer- 
romagnetism according to the S = I and S = | high-spin 
M cations, respectively 16 ' 17 ' 18 For these delafossites, the 
complex magnetism resulting from the frustrated nature 
of the MO2 network is indeed thought to be responsible 
for the magnetic field induced electric polarization. 

Such multiferroic behavior is in marked contrast to the 
metal-type behavior reported for the Mg 2+ substituted 
rhodate CuRhi_ x Mg x 02. The 4d 5 /4d 6 electronic config- 
urations of low-spin Rh 4+ /Rh 3+ in this series provides 
a unique opportunity to study the electronic ground- 
state of a delafossite without extra contributions gen- 
erated by the large magnetic moments of Cr 4+ /Cr 3+ 
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(S = 1/5 = §) or Fe 4+ /Fc 3 + (S = 2/S = §) cations. 
However, literature data for CuRh 1 _ x Mg x 2 show some 
variation. This is best illustrated by the different values 
reported for the room temperature Seebeck coefficient S 
of CuRh .9oMgo.io0 2 , which are 5 30 ok = 130/xVK -1 , 
and S300K = 70^VK _1 as given in Refs. H and 0, 
respectively.— Among the possible reasons explaining 
this discrepancy is the uncertainty of the Mg for Rh sub- 
stitution controlling the concentration of Rh 4+ holes in 
the Rh 3+ matrix. In order to check for the transport 
mechanism and for the substitution effectiveness, elec- 
tronic structure calculations have been performed. In 
the present paper, these results are discussed in com- 
parison to the experimental data. The latter are ob- 
tained on polycrystalline samples, for which the delafos- 
site structure was verified by a combined X-ray and elec- 
tron diffraction study. Taken together with the cation 
analysis, these characterizations also demonstrate the ex- 
istence of a Mg 2+ solubility limit much larger than in 
CuCri_ x Mg x 02. The electrical resistivity and Seebeck 
coefficient measurements show that the substitution of 
Mg for Rh in CURI1O2 allows to progressively induce an 
insulator to metal transition accompanied by a gradual 
decrease of the thcrmopower. 

Contrasting previous studies, these results arc found to 
be compatible with electronic structure calculations with 
a lack of significant contribution of copper, the most im- 
portant contribution at the Fermi level coming from the 
Ad orbitals of the Rh cations. The lack of significant mag- 
netic contribution, probed by magnetic susceptibility, re- 
sponsible for the low temperature increase in CuCr02 
and CuFe02, allows for the study of the metallic state 
from 2.5 K to 1000 K. An unusual T 2 regime is evidenced 
over a broad temperature range, which, when combined 
with the T dependence of the Seebeck coefficient, leads 
to remarkably T independent power factors, PF = — . 



II. METHODOLOGY 

A. Electronic structure calculations: theoretical 
method 

The calculations are based on density-functional 
theory and the generalized gradient approximation 
(GGA)3 They were performed using the scalar- 
relativistic implementation of the augme nted spherical 
wave (ASW) method (see Refs. HHHl and references 
therein). In the ASW method, the wave function is 
expanded in atom-centered augmented spherical waves, 
which arc Hankel functions and numerical solutions of 
Schrodingcr's equation, respectively, outside and inside 
the so-called augmentation spheres. In order to optimize 
the basis set, additional augmented spherical waves were 
placed at carefully selected interstitial sites. The choice 
of these sites as well as the augmentation radii were au- 
tomatically determined using the sphere-geometry opti- 
mization algorithm^ Self-consistency was achieved by a 



highly efficient algorithm for convergence acceleration^ 
The Brillouin zone integrations for the self-consistent 
field calculations were performed using the linear tetra- 
hedron method with up to 1156 k-points within the ir- 
reducible wedge of the rhombohedral Brillouin zone ) 23 i 26 
whereas the calculation of the density of states and the 
thcrmopower is based on 19871 k-points. 

In the present work, a new full-potential version of the 
ASW method was employed^ In this version, the elec- 
tron density and related quantities arc given by spherical- 
harmonics expansions inside the muffin-tin spheres. In 
the remaining interstitial region, a representation in 
terms of atom-centered Hankel functions is used.— How- 
ever, in contrast to previous related implementations, we 
here get away without needing a so-called multiple-K ba- 
sis set, which fact allows to investigate rather large sys- 
tems with a minimal effort. 



B. Ceramic samples preparation and 
characterization 

The polycrystalline samples of the CuRhi_ x Mg x 02 se- 
ries have been prepared by solid state reaction in air. 
Bars of typical size 2 x 2 x 10 mm were prepared by 
mixing stoichiometric amounts of the CU2O, RI12O3 and 
MgO precursors, which were then pressed. The electron 
diffraction study was carried out with a JEOL 201 OCX 
transmission electron microscope. The sample prepara- 
tion is made by crushing in butanol some bar pieces, the 
corresponding microcrystals being afterwards deposited 
on Ni grids. The purity of the obtained black bars was 
checked by X-ray powder diffraction using a Panalytical 
X-pert Pro and a Brucker diffractometers. The data were 
analyzed by using the Fullprof suited 

As two different calcination temperatures have been 
used in the literature, two attempts were made at ei- 
ther 930°C or 1050°C for a duration of twelve hours. 
As shown in Fig. Q}t for CuRho.9Mgo.1O2 (x = 0.1), one 
small extra peak identified as CuO appears (see the ver- 
tical arrow in Fig. QJi) . The weight of the latter depends 
on the synthesis temperature. Indeed, when increasing 
it from 930°C to 1050°C, the weight of this peak de- 
creases. Furthermore, as these patterns were recorded in 
the same conditions, (amount of powder and acquisition 
time) , this comparison shows the poor crystallinity of the 
930°C prepared sample. This result agrees with the less 
dispersed cation contents of the 1050°C prepared sam- 
ple as probed by energy dispersive X-ray spectroscopy 
(EDS) analysis performed within the transmission elec- 
tron microscope. Such conclusions are consistent with 
the thermodynamics of the Cu-Rh-0 ternary diagram, 
which showed that CuRh02 synthesized below 985°C in 
an oxygen pressure of 0.1 MPa is unstable^ Thus, the 
sample scries corresponding to x = 0, 0.01, 0.04, 0.10, 
0.15, and 0.30 in CuRhi_ x Mg x 02 has been calcined in 
air at 1050°C for a duration of twelve hours. 

The patterns of the CuRhi_ x Mg x 02 samples are all 
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The low temperature (T < 320 K) electrical resistiv- 
ity (p) and Seebeck coefficient (S) were measured by 
using a Quantum Design cryostat. The four-probe and 
the steady-state techniques were used for the former and 
the latter, respectively, indium electrical contacts having 
been deposited with ultrasons. For the high temperature 
p and S measurements (T > 300 K) a Ulvac-Zem 3 sys- 
tem was used. The magnetic susceptibility was measured 
by using a dc SQUID magnetometer (zero field-cooling, 
p Q H = 0.3 T). 
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FIG. 1: (Color online) Top: X-ray patterns of two 
CuRho.goMgo.10O2 samples synthesized at 930°C and 1050°C, 
respectively, with indexation in the R3m space group. Bot- 
tom: X-ray diffraction patterns at room temperature; the 
vertical bars correspond to the locations of the diffraction 
peaks in the R3m space group (with a — 3.0741(1) A, and 
c = 17.0952(3) A). The missing data in the 44-44.8° region 
(29) have been intentionally eliminated due to a small peak 
coming from the sample holder. 

refined in the R3m space group, usually reported for de- 
lafossite compounds at room temperature, as shown in 
Fig. [U A small decrease of both a and c lattice parame- 
ters is observed when x increases up to 15% (correspond- 
ing to a decrease of ~ 0.6% of the cell volume). From 
x = 0.10 on, one small extra peak appears that can be 
identified as CuO. 

In order to test the maximum content of Mg 2+ substi- 
tuted for Rh cations, EDS analyses coupled to electron 
diffraction have been made first for the compound with 
the highest experimental Mg doping (x = 0.3). This ana- 
lytical study demonstrates that a maximum of 12% Mg 2+ 
can be substituted at the Rh site in CuRhC>2, while the 
Cu2Mg03 oxide is detected as a secondary phase. This 
is consistent with the observation from X-ray diffraction 
of this impurity for x > 0.10. For that reason, our mea- 
surements of transport properties are restricted to com- 
pounds corresponding to x < 0.10, the x = 0.10 compo- 
sition already containing a very low amount of CuO as 
impurity. We emphasize that this structural study can- 
not be reconciled with previous data by Shibasaki et al4 
given for compounds nominally containing up to 20% of 
Mg 2+ . 



III. RESULTS 

A. Calculations: Important role of the Rh 4c? 
orbitals 

The calculations were based on the crystal structure 
data by Oswald et aZ-ri who determined the lattice con- 
stants as a — 3.08 A and c = 17.09 A. However, these 
authors did not measure the internal oxygen parameter. 
For this reason, we performed an energy minimization, 
leading to a value of zq = 0.10717, which was used in all 
subsequent calculations. 

The electronic bands along selected high-symmetry 
lines of the first Brillouin zone of the hexagonal lattice, 
Fig. [2j are displayed in Fig. [3] 




FIG. 2: First Brillouin zone of the hexagonal lattice. 



The corresponding partial densities of states (DOS) are 
shown in Fig. [4] While the lower part of the spectrum is 
dominated by O 2p states, the transition metal d states 
lead to rather sharp peaks in the interval from —4 to 
+4eV. In particular, the ti g and e g manifolds of the Rh 
Ad states as resulting from the octahedral coordination 
are recognized. This representation of the partial DOS 
used a local rotated coordinate system with the Cartesian 
axes pointing towards the oxygen atoms, c-type overlap 
of the O 2p states with the Rh Ad e g orbitals leads to the 
contribution of the latter between —7 and — 6eV. In con- 
trast, due to the much weaker 7r-type overlap of the O 2p 
states with the t2 g orbitals, these states give rise to sharp 
peaks in the interval from — 2.2eV to the valence band 
maximum. The ti g manifold is separated by an optical 
band gap of ps 0.75 cV from the empty e g states and thus 
leads to a Rh d 6 state. The Cu 3d states are essentially 
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limited to the interval from — 4eV to the valence band 
maximum and thus Cu can be assigned a monovalent 
d 10 configuration in close analogy with the experimental 
findings. In passing, we mention the finite dispersion of 
the electronic bands parallel to T-A, which points to a 
considerable three-dimensionality arising from the cou- 
pling between the layers. Yet, this overall behavior of 
the dispersion perpendicular to the a — b plane is con- 
trasted by the barely noticeable dispersion particularly 
along the line K-H. This has been also observed for other 
dclafossitc materials . 10 ' 12 

Further insight about the electronic properties of 
CuRh02 can be gained by analyzing the real and imagi- 
nary parts of the dielectric function as calculated within 
linear-response (see Ref. for more details). As is ev- 
ident from Fig. the asymmetry between the in-plane 
and out-of-plane directions is not reflected in the absorp- 
tion gap following from the imaginary parts of the di- 
electric function. In fact the gap is very close to 0.75 eV 
in all three directions, which is most likely to exceed the 
Hund's rule coupling. Therefore, the low-spin Ad 6 config- 
uration of Rh 3+ is expected to be the ground state. This 
is consistent with earlier findings by Singh for CUC0O2, 
where the Co ions adopt the low-spin 3d 6 configuration. 8 

Finally, we have calculated the thermopower using the 
framework of Boltzmann theory.— The transport prop- 
erties are expressed in terms of the Onsager transport 
coefficients, 
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(i) 

where ( — ^ is the negative derivative of the Fermi 
function and 

Saa< (E) = -j- E < v L^nS(E - e kn ) (2) 



k n 



denotes the so-called transport distribution i 33 ' 34 Here, fi c 
is the volume of the unit cell, v£ n a Cartesian component 
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FIG. 3: (Color online) Electronic bands of CuRh02- 
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FIG. 4: (Color online) Partial densities of states (DOS) of 
CuRh02- Selection of the Rh 4c? orbitals is relative to the 
local rotated reference frame, see text. 




FIG. 5: (Color online) Dielectric function of CuRh02 



of the group velocity of the ri'th band, and T kn is the 
relaxation time. While the electrical conductivity and 
the thermal conductivity (at zero electric field) are given 
directly by JT]) for n = and n = 2, respectively, the 
thermopower is calculated from the matrix equation 



S 



A.V 



1 

eT 



L (o) 
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XX' 



where e is the (negative) electronic charge. Following 
standard practice we assume that the relaxation time 
does not depend on k and band index, in which case r 
cancels from the Sccbeck coefficient. The implementation 
of the above formulation was done along similar lines as 
those previously proposed by Scheidcmantcl et al. as well 
as by Madsen and Singh. 33 ' 34 The implementation was 
tested against the recent results by Singh for CUC0O2 
and YCu02 and very good agreement was founds 

The xx-components of the thermopower as calculated 
for different doping levels are displayed in Fig. the 
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FIG. 6: (Color online) Thermopower S xx of CuRh02 for dif- 
ferent hole doping levels. 



calculated zz-components are about 10-20% larger. Ac- 
cording to these findings, the thermopower strongly de- 
creases with increased hole doping. In addition, it shows 
an almost linear dependence on temperature especially in 
the intermediate-temperature range with the downturn 
at low doping and high temperatures reflecting excita- 
tions across the optical band gap. Worth mentioning arc 
the rather high values for the lower doping levels down 
to ?s 100 K and the pronounced drop below this temper- 
ature. 

In passing, we mention that apart from systematically 
slightly smaller values our results are in perfect agree- 
ment with the calculations of Usui et al.^ who likewise 
used the Boltzmann equation approach and who in turn 
obtained almost perfect agreement with the experimen- 
tal data by Kuriyama et al.£ However, we recall from 
the above mentioned previous comparative tests to the 
results by Singh that the thermopower is remarkably sen- 
sitive to details of the crystal structure. Since Usui et al. 
used slightly different lattice constants in their calcula- 
tions this might explain the systematic deviations in the 
calculated thermopower. 



B. Electrical resistivity: A metal-insulator 
transition 



p(T) curves for the CuRhi_ x Mg x 02 series as given in 
Fig. [7] reveal the modification of the electronic ground- 
state induced by the substitution. At 300 K, the val- 
ues decrease by a factor of ~ 300 as x increases from 
x = 0.00 to x = 0.10. This drop is even more pronounced 
at lower temperatures as the p(T) curve exhibits a local- 
ized behavior in CuRh02 with Piook = 3 10 4 f2 cm at 
100 K. In contrast, metal-like behavior is observed for 
CuRho.90Mgo.10O2 with piooK = 1-2 10 _3 fi cm. It must 
be emphasized that the change of electronic state induced 
by Mg 2+ is progressive. The CuRho.99Mgo.01O2 com- 
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FIG. 7: (Color online) Temperature dependence of the resis- 
tivity of CuRhi_ x Mg x 2 for x = 0, x = 0.01, x = 0.04, and 
x = 0.10. 



pound still exhibits a localized behavior but with a resis- 
tivity decreased by three orders of magnitude at 100 K as 
compared to CURI1O2, whereas for CuRho.96Mgo.04O2, 
p remains almost T independent from 100 K to 1000 K 
(p increasing by 5% only in this temperature range). 
The p(T) curve shows a re-entrant behavior only below 
~ 100 K, reaching a maximum value of p = 8 10 _1 Q cm 
at 5K. In fact, a closer inspection of the curves re- 
veals that they all go through a minimum value at a 
characteristic temperature T m i n separating a -Jp < 
regime below T m - m from a > regime for T > T m ; n 
(Fig. [7J. As shown in Table [H the T m ; n value decreases 
from T min = 800 K for CuRh0 2 to T min = 38 K for 
CuRho.90Mgo.10O2. 

For the pristine compound, the localizing behavior 
observed below T m i n (800 K) is consistent with the ex- 
istence of a rather small gap at the Fermi level ob- 
tained in the section IIII Al Even though this tem- 
perature dependence of the resistivity is similar to 
the curve measured for CuCr02, it is in fact closer 
to the one of CuCro.9gMgo.oi02»^ However, the data 
point to a different transport mechanism. Indeed, first, 
no significant magneto-resistance was observed in con- 
trast to all CuCi'i_ x Mg x 02 samples showing a magneto- 
resistance as high as —10% at 5K in 7T. Second, 
neither the Arrhenius law p oc e~ T °/ T nor the pola- 
ronic model p cx Te~ T °/ T , nor the variable range hop- 



x 


T min (K) 


p(T , min)(mf2 cm) 


0.00 


800 


102 


0.01 


699 


45.56 


0.04 


320 


10.74 


0.10 


38 


1.22 



TABLE I: T m i n and doping-dependence of the resistivity at 
T m i n for several compositions of CuRhi_ x Mg x 02. 
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ping model p oc e"' T °/ T ' with a = 1/2, 1/3 or 1/4, 
which are broadly used in conventional three dimen- 
sional transition metal perovskites as Lai_ x Sr x Co03,— 
and which were successfully applied to the two di- 
mensional chromium based delafossites CuCr02 and 
CuCro.98Mgo.02O2, respectively— can convincingly fit 
the p(T) data of CuRh0 2 . 

This conclusion about the transport mechanism in 
CuRh02 is confirmed by the analysis of the p(T) curves 
found for the Mg-substituted CuRh02 compounds. 
First, the attempts to fit the curves by a polaronic model 
fail as for CuRh02. Second, a low temperature Fermi liq- 
uid behavior is found for CuRho.90Mgo.10O2. Third, the 
metal-like regions of the curves can be adjusted to T 2 
dependences that are obeyed over a wide temperature 
range. As shown in Fig. [8] for CuRho.96Mgo.04O2, the 




°0 20 40 60 80 100 

2 4 2 

T (10 K ) 

FIG. 8: (Color online) Fermi liquid-like behavior of 
CuRhi-xMgxCb for the dopings x = 0.04 and x = 0.1. For 
CuRho.g6Mgo.04O2, a constant contribution p cc = 8mf2 cm 
has been subtracted for clarity. The thin full lines represent 
the Fermi liquid fits, while the thick full lines are guides for 
the eyes only. 

p oc T 2 regime holds from ~ 300 K to ~ 850 K. Express- 
ing p(T) as p(T) = po + AT 2 as in a Fermi liquid we 



Compound 


A(Q m K-' 2 ) 


CuRho.g6Mgo.04O2 
CuRho.goMgo.10O2 
PdCo0 2 (from Ref. Ft) 
Thin films of V2O3 (from Ref. 36) 
Lai_ x Sr x Ti0 3 (from Ref. 37) 


7.7 10 _il 
3 10" 11 

4.8 10~ 13 
2 10" 9 

2-3 10" 11 



TABLE II: Fermi liquid transport parameter for several ox- 
ides. 



obtain the transport parameter A. As shown in Table HT1 
the values of A for these rhodates are found to be lo- 
cated in an intermediate range, namely they arc larger 
than in the PdCo02 delafossites^ but smaller than in 
thin films of V2O312& In a fashion similar to the behavior 
observed in the titanates Lai- x Sr x Ti03,— A increases 
with decreasing concentration of charge carriers, while, 
in contrast to all the above systems, the T 2 -behavior may 
be observed for temperatures up to 1000 K. Remarkably 
the widely observed phonon-dominated p oc T behavior 
is not taking over, even at such high temperatures. Such 
a result enlights the unusual transport behavior of these 
rhodates. The presence of the Mg 2+ scattering centers, 
up to 12%, corresponds to concentration well below the 
percolation threshold, and no band purely based on Mg 
orbitals is expected to form. Therefore, on their own, 
such low Mg concentrations should not affect the trans- 
port in the T range where this p oc T 2 regime is observed 
but by hole doping the Rh-based Ad bands. Still, inho- 
mogeneous distribution of the Mg ions on the Rh sites 
might be responsible of the kink observed in Fig. [8] for 
CuRho.90Mgo.10O2. 

Finally, the resistivity values for CuRhi_ x Mg x 02 are 
rather comparable to those reported in Ref.[j (T < 300 K) 
or in Ref. y (T > 400 K). The decrease of p induced by 
the Mg 2+ substitution in the present samples strongly 
suggest that "hole" charge carriers are created according 
to the formula CuRhit 2x Rh x + Mg 2 +0 2 . It must also be 
added that since no other CuM0 2 delafossites exhibits 
such metal-like behavior down to very low T, the role of 
the Cu channel to the electronic transport can hardly be 
invoked as is also confirmed by the electronic structure 
calculations. 



C. Thermoelectric power and power factors 

Although the p values of our CuRhi_ x Mg x 02 series are 
comparable to those already reported^ the Seebeck coef- 
ficients values for CuRh02 (Table ITTT1 and Fig. [9]) appear 
to be different. When compared to the study in which 
the samples were calcined at lower T ^ the present S val- 
ues are found to be always much larger, as documented 
in Table [TTT1 Besides, our data for the substituted com- 
pounds showing an S decrease as x increases cannot be 
reconciled with the x-independent S300K for all substi- 
tuted compounds of Ref. |4|. 

For the most metallic sample, the S values are al- 



x 


T(K) 


S(pV KT 1 ) 





330 


280 


(from Ref. 4) 


300 


130 


0.10 


450 


120 


0.10 (from Ref. 3) 


450 


165 


0.05-0.20 (from Ref. £j 


450 


65 



TABLE III: Typical values of the thermopower of 
CuRhi_ x Mg x 2 . 
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FIG. 9: (Color online) Temperature dependence of the ther- 
mopower of CuRhi- x Mgx0 2 for x = 0, x = 0.01, x = 0.04, 
and x = 0.10. Lines are guide for the eyes, only. 



ways increasing with T as shown for CuRh .goMgo .10O2 
in Fig. [91 In contrast, the localizing behavior of CuRh02 
is reflected by the S(T) curves showing an upturn to- 
wards high S values as T decreases below ~ 450 K. For 
all x > 0.01, or for T > 450 K (for x = 0.00 and x = 
0.01) the S(T) curves exhibit an almost S oc T regime as 
shown by the lines drawn on the S(T) curves in Fig. [9] 
Such a T dependence is characteristic of a metallic be- 
havior. As previously reported for CuRho.90Mgo.10O2, 
this behavior leads to rather large positive values, larger 
than 300 /iV K" 1 at 1000K for x = 0.00 and x = 0.01. 
The decrease of S as x increases together with the pos- 
itive sign of S is consistent with an increase of the hole 
(Rh 4+ ) fraction induced by charge compensation created 
by the Mg 2+ for Rh 3+ substitution. 

Remarkably, all these trends are in good qualitative 
agreement with the theoretical findings presented in Sec. 
IIII Al Yet, the experimentally obtained values at 1000 K 
are somewhat smaller than the theoretically predicted 
ones. One is therefore tempted to analyze the above ex- 
perimental data within a completely different approach, 
namely, the temperature independent correlation func- 
tions ratio approximation (TICR)i^ In this approach, 
the thermopowcr arising from the Kubo formalism 



S(T) 



1 {jEjn) - fJ-{jnjn) 



eT 



(jnjn) 



(4) 



is approximated by assuming that the energy current- 
particle current correlation function and the particle cur- 
rent autocorrelation function share the same temperature 
dependence. Thus their ratio should result in a hyper- 
bolic offset of strength Eq of the thermopower as 



S(T) = -L(E -»(T)). 



(5) 



At high temperatures the temperature dependence of the 
thermopower is governed by the one of the chemical po- 



n(T) = / def(E - fi)p(E) 



(0) 



Here, f(E) again denotes the Fermi function and p(E) is 
the density of states. In the present context, the latter is 
taken from the electronic structure calculations (Fig. 0|. 
Yet, within the present model there is still room for im- 
provement especially at high temperatures. This goes 
mainly along two different directions. One way to im- 
prove on the TICR as covered by I©/© would be to 
use a modified density of states resulting from the GGA 
result by a rigid energetical upshift of the conduction 
bands. However, for reasonable values of the latter, it 
turned out that the thermopower shows little sensitivity 
to the actual optical band gap. From this, we can fur- 
thermore conclude that the behavior of the thermopower 
is dominated by the holes in the valence bands. 

Another direction is provided by the observation that, 
in many cases, the thermopowcr tends to loose its tem- 
perature dependence above room temperature, see e. g. 
Ref. [39|. With this motivation, the TICR expression (J2|) 
was extended by adding a temperature independent con- 
tribution Sq to S(T). Such a contribution is often at- 
tributed to localized degenerate states. Here they may 
derive from the pockets on the Fermi surface centered 
around the H and A points that are characterized by 
small Fermi velocities (see Fig. [3]). As a consequence, 
within the extended TICR the thermopower is deter- 
mined by two parameters, namely, Eq and So, which, 
together with the GGA DOS, can be used to fit the ex- 
perimental data. The result is shown in Fig. [10] for all 
three doped samples. Regarding Eq, the values range 
from 40 mcV to 90 meV, thereby being quite similar 
to the ones reported for electron-doped manganites. In 
these materials a large density of states at the bottom 
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FIG. 10: (Color online) Comparison of the different theo- 
retical models, including the temperature independent con- 
tribution So, (lines) with our experimental values (sym- 
bols) of the temperature dependence of the thermopower of 
CuRhi_ x Mg x 2 for x = 0.01, x = 0.04, and x = 0.10 
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of the band results in low degeneracy temperatures and 
large negative thermopower^ For the title compounds, 
the density of states at the top of the valence band is 
large, resulting in the large positive thermopower, to- 
gether with low degeneracy temperature. The So values 
are quite large (AO^VR- 1 < S < lOOpVK- 1 ) and ap- 
pear to be comparable to the ones reported for various 
layered colbatatcs, in which the here considered RhC>2 
layers are replaced by isostructural C0O2 layers^ Addi- 
tionally, the fact that the thermopower data for all three 
doped samples could be reproduced using the same DOS 
for both the extended TICR model and the GGA calcu- 
lations is quite remarkable. This gives a strong support 
to use the rigid band model. 

In order to check the thermoelectric performance, we 
now address the power factor PF (PF = — ). Combin- 
ing their T dependences, S oc T and p oc T 2 at suf- 
ficiently high temperature, this leads to remarkable T- 
indcpcndent values of the PF as shown in Fig. [TT] The 
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FIG. 11: (Color online) Temperature dependence of the power 
factor of CuRhi_ x Mg x 2 for x = 0, x = 0.01, x = 0.04, and 
x = 0.10. Lines are guide for the eyes, only. 

best values are observed for the most substituted sam- 
ples exhibiting a value of the PF = 7 10 _4 WK _2 m" 1 . 
As shown in Fig. [TTJ and in good qualitative agreement 



with the calculations presented in Sec. IIII A[ the values 
are found to increase as x increases in CuRhi_ x Mg x 02 
showing that the induced relative p decrease is more than 
compensating the S decrease. 



IV. CONCLUSION 

In summary, Mg-doped CURI1O2 has been investigated 
by means of electronic structure calculations, structural 
characterization, and transport measurements. The elec- 
tronic structure calculations clearly indicate that the 
transport is dominated by the Rh 4c? bands. Struc- 
tural data demonstrate that the solubility limit of Mg 
in CuRh02 is as high as 12%, provided the samples are 
prepared at temperatures above ~ 1000°C, in which case 
Rh is indeed substituted by Mg. This substitution re- 
sults in a peculiar hole doping of the rather narrow Rh 
Ad bands, shown by a T 2 dependence of the resistiv- 
ity, a T dependence of the thermopower, and a quite 
large nearly T independent power factor, up to temper- 
atures as high as 1000 K. Regarding the thermopower, 
good qualitative agreement between the theoretical pre- 
diction arising from the GGA+Boltzmann approach and 
experimental data is obtained. Yet, an additional purely 
cntropic contribution needs to be invoked when treating 
the GGA results in the TICR framework for quantita- 
tive explanation, in a similar fashion as what was found 
in cobaltites with isostructural C0O2 layers.— Thermal 
conductivity measurements performed on dense samples 
would be necessary to measure the figure of merit. Work 
along this line is in progress. 
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